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Abstract 

We present the first clear observation of the doubly Cabibbo suppressed decay D + — > 
K~K + K + and the first observation of the singly Cabibbo suppressed decay Df — > 
K~K + K + . These signals have been obtained by analyzing the high statistics sample 
of photoproduced charm particles of the FOCUS (E831) experiment at Fermilab. 
We measure the following relative branching ratios: 

T (D + -► K-K+K+) /r (D+ -»• K--K+-K+) = (9.49 ± 2.17 ± 0.22) x 1(T 4 

and 

T (Df -► K-K+K+) /r (D+ -► K-K+7T+) = (8.95 ± 2.12 t\ H) x 10" 3 , 
where the first error is statistical and the second is systematic. 



1 Introduction 



Doubly Cabibbo suppressed (DCS) charm decays are expected to occur with a 
rate which is roughly a factor tan 4 9c ~ 2.5 x 10 -3 smaller than the correspond- 
ing Cabibbo favored (CF) modes. This is the main reason our present knowl- 
edge of these decays is rather poor and limited to very few decay modes. Only 
four DCS decays have been observed, D + — > K + -n~-n + , D° — > K + n~, K + it~it 



* See http://www-focus.fnal.gov/authors.html for additional author information. 
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and K + h~'k + 'k~ .Q The interpretation of the D° modes is complicated by pos- 
sible contributions from D° — D° mixing [4-7], making the D + — > K + tt~tt + 
decay the only pure DCS decay previously studied. 

In this paper, we report the first clear observation of the DCS decay D + — > 
K~K + K + , together with the first observation of the singly Cabibbo sup- 
pressed (SCS) decay of Df into the same final state. Throughout this paper, 
the charge conjugate is implied when a decay mode of a specific charge is 
stated. 



It is interesting to note that in contrast to the four modes previously men- 
tioned, the DCS decay D + — > K~K + K + cannot result from a simple spec- 
tator process, but presumably requires the intervention of strong resonances 
that simultaneously couple to the tttt and KK channels. It could also pro- 
ceed through annihilation but from studies of Df — > we expect this 
contribution to be small [8]. 

The results presented in this paper have been obtained using the high statis- 
tics charm sample of the FOCUS experiment at Fermilab. FOCUS is a charm 
photoproduction experiment which took data during the 1996/1997 fixed tar- 
get run at Fermilab. The FOCUS detector is a large aperture, fixed-target 
spectrometer with excellent vertexing and particle identification. A photon 
beam is derived from the bremsstrahlung of secondary electrons and positrons 
with an rj 300 GeV endpoint energy produced from the 800 GeV/c Teva- 
tron proton beam. The photon beam interacts in a segmented BeO target. 
The charged particles which emerge from the target are tracked by two sys- 
tems of silicon microvertex detectors. The upstream system, consisting of 4 
planes (two views in 2 stations), is interleaved with the experimental target, 
while the other system lies downstream of the target and consists of twelve 
planes of microstrips arranged in three views. These detectors provide high 
resolution separation of primary (production) and secondary (decay) vertices 
with an average proper time resolution of ~ 30 fs for 2-track vertices. The 
momentum of a charged particle is determined by measuring its deflections 
in two analysis magnets of opposite polarity with five stations of multiwire 
proportional chambers. Three multicell threshold Cerenkov counters are used 
to discriminate between electrons, pions, kaons, and protons. 



1 Evidence for the DCS decay D + — > K~ K + K + was previously reported by two 
experiments [1,2], but their results were superseded [14] by the much more stringent 
upper limits coming from the higher statistic experiment E687 [3]. 



3 



2 Signals and selection criteria 



The final states are selected using a candidate driven vertex algorithm. The 
basic idea of this algorithm is to use a charm candidate decay vertex as a seed 
to find the primary vertex. In our particular case a decay vertex is formed from 
three reconstructed charged tracks and the momentum vector of the resultant 
D candidate is used to intersect other reconstructed tracks and search for a 
suitable production vertex. The confidence levels of both vertices are required 
to be greater than 1%. We measure I the separation of the two vertices and its 
associated error at. The quantity tjot is the significance of detachment of the 
secondary and primary vertices. Cuts on tjat are used to extract the D signals 
from non-charm background and to improve the signal to background ratio. 
Two other measures of vertex isolation are used: a primary vertex isolation 
and a secondary vertex isolation. The primary vertex isolation cut requires 
that the confidence level for one of the tracks assigned to the decay vertex to 
be included in the primary vertex be less than a certain threshold value. The 
secondary vertex isolation cut requires that the maximum confidence level for 
all tracks not assigned to any vertex to form a vertex with the D candidate be 
less than a certain threshold value. The main difference in the selection criteria 
between different decay modes lies in the particle identification cuts applied to 
the decay products. To minimize the systematic errors we use identical vertex 
cuts both on the signal and normalizing modes. 

In the D + — > K~K + K + analysis we require £/a e > 8. The primary and 
secondary vertex isolation must be less than 0.1%. The D momentum must be 
in the range 25 GeV/c to 250 GeV/c and the primary vertex must be formed 
with at least two reconstructed tracks in addition to the seed track. We require 
that the decay vertex occur outside of the target material. For each charged 
track the Cerenkov algorithm computes four likelihoods from the observed 
firing response of all the cells that lie inside the track's Cerenkov cone for every 
counter [9] . The product of all firing probabilities for all cells within the three 
Cerenkov cones produces a x 2 -like variable Wi = — 2 ln(Likelihood), where % 
ranges over electron, pion, kaon and proton hypotheses. We require observed 
Cerenkov light pattern for the kaon hypothesis is favored over that for the pion 
hypothesis by more than a factor of exp(0.5) by requiring W w — W K > 1.0. We 
also apply a kaon consistency cut, which requires that no particle hypothesis 
is favored over the kaon hypothesis with a AW = Wk — W m [ n exceeding 3.5. 
To further reduce the background due to poorly reconstructed candidates, we 
require that the proper time resolution of the candidates, defined as cre/^jc), 
be less than 150 fs. 

The resulting D + signal is shown in Fig. 1(a). We obtain a Gaussian yield 
of 65.5 ± 15.0 D + — >■ K~K + K + events over a linear background. The mass 
value returned by the fit is 1869 ± 1 MeV/c 2 ; the r.m.s. of the Gaussian fit is 
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5.2 ± 1.2 MeV/c 2 in agreement with Monte Carlo simulations. The two broad 
structures around 1985 MeV/c 2 and 2085 MeV/c 2 are due to D + and D+ 
decays into K~K + tt + where the n + is misidentified as a K + . 

In the Df — >■ K~K + K + analysis we have to use stronger Cerenkov cuts 
to extract the signal which otherwise would be completely hidden by the 
K~ K + -k + mis-identification peaks. We require W n — Wk > 4.5 for all three 
kaon candidates. All the other cuts are the same as for the D + — > K~K + K + 
decay. 

Fig. 1(b) shows the invariant mass plot where both D + and peaks are now 
evident. In the fit the mass and width are fixed to the values found in the 
Monte Carlo. This is done to reduce the effects of any residual fluctuation of 
the D + — > K~K + ir + reflection, which would induce a shift of the peak toward 
higher masses. We obtain a yield of 31.4 ± 7.4 — > K~K + K + events over 
a linear background. 

For D + — > K~K + K + we measure the branching ratio relative to D + — > 
K-ir + iT + , while for D+ -> K~K + K + that relative to D+ -> K~K + n + . We 
obtain: 

T (D+ -> K~K + K + ) /r (D+ -> fsr-7r + 7r+) = (9.49 ± 2.17) x 10" 4 

T (D+ -> K-K+K+) /r (D+ -> X+K-7T+) = (8.95 ± 2.12) x lO" 3 . 

The cuts on the normalization modes are identical whenever possible to those 
used for the selection of the corresponding 3K signal. In addition, to remove 
contamination from the — > K~K + tt + normalization mode due to Cerenkov 
misidentified D + — > K~h + 'k + events, we employ an anti-reflection cut to 
reject candidates which, when reconstructed as K~it + -n + , lie within 2 sigma 
of the D + nominal mass. The normalization signals are shown in Fig. 1(c) and 
Fig. 1(d) and consist of 62911 ± 263 and 3844 ± 66 events respectively. 

In all our simulations we always used the proper resonant substructure for the 
two normalization modes [10] [11], which would otherwise produce important 
systematic deviations of the results. 



3 Systematic Errors 

We performed a detailed investigation of any source of systematics which could 
impact our branching ratio measurements. We first studied the stability of the 
results by varying the cuts over a wide range of values. Our results are stable 
in their evolution on the most critical cuts: £/&£, W n — Wk and primary and 
secondary vertex isolation. 
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1.7 1.8 1.9 2 2.1 1.7 1.8 1.9 2 2.1 

KK + K + Mass GeV/c 2 K'K + K + Mass GeV/c 2 



(a) (b) 



rT- 1400 




1.7 1.8 1.9 2 2.1 1.7 1.8 1.9 2 2.1 

K'jiVMass GeV/c 2 K"KV~Mass GeV/c 2 



(c) (d) 

Fig. 1. Invariant mass distributions for D+ -> K~K + K + {a), D+ -> i^K+K+tV), 
£>+ -► K-7T+7T+(c) and L>+ -> K-K+7T+(d). 

We then split the samples using variables which can probe different kinemati- 
cal regions, such as low and high momentum range, or different experimental 
conditions, such as early and late runs, which have different target configura- 
tions. In doing this we can check our results together with our Monte Carlo 
simulation over a variety of different conditions. We quantify a "split sample 
systematic error" by examining consistency among these statistically indepen- 
dent splits of our data. If the consistency \ 2 turns out to be smaller than f , 
this error is taken to be zero. Otherwise we scale all the errors up to bring the 
X 2 back to f . The split sample systematic error is then defined as the difference 
in quadrature between the scaled error of the weighted average of the subsam- 
ple estimates and the statistical error of the total data set. This procedure is 
similar to the S'-factor method used by the Particle Data Group [14]. 

We have split our sample by high and low D-momentum, D and D, and early 
and late run periods. Splits have been done in one variable at a time because 
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of our limited statistics. 



The measured branching ratios for the three pairs of disjoint samples are shown 
in Fig. 2. We find only one contribution to the systematic uncertainty, namely 
the run-period split sample for the decay which gives a contribution to 
the branching ratio systematics of 2.23xl0~ 3 . 



x 10 
0.14 



0.02 - 
0.0175 - 
0.015 - 
0.0125 - 



Split Samples 



Split Samples 



(a) 



(b) 



Fig. 2. Split sample results for -D(a) and D s (b) relative branching ratios. Three pairs 
of disjoint samples are considered: high and low momenta on the left, late and early 
runs in the center, D and D on the right. The lines show the joint sample and the 
la error bars. 

In computing the branching ratios we have used the efficiency of a pure phase- 
space decay. This choice was motivated by the relatively flat distribution of the 
events over the Dalitz domains as shown in Fig. 3. To better investigate the 




1.4 1.6 1.8 2 1 1.2 1.4 1.6 1.8 2 2.2 

m(KK) 2 (GeV/c 2 ) 2 m(KK) 2 (GeV/c 2 ) 2 



(a) (b) 



Fig. 3. Dalitz plot for D + (&) and for Df(h). Only events which lie within 2a of the 
respective nominal masses are plotted. 
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implications of this assumption we have computed the reconstruction efficien- 
cies for two particularly representative (pK + decay and a /o(980)AT + 
decay. Table 1 shows the calculated efficiencies with respect to those for pure 
phase-space decays. Given the non-negligible variation of the efficiency val- 





e(D+) 


e(Df) 


Phase-Space 


1 


1 


cj)K + 


0.927 ±0.015 


0.948 ± 0.015 


/ (980)^+ 


1.028 ±0.014 


1.086 ± 0.014 



Table 1 

Reconstruction efficiencies, e, for different decay dynamics into the same K~ K + K + 
final state for D + and Df. 

ues, we considered the following two cases in order to assess the systematic 
uncertainty: the decay proceeds through the maximum estimated amount of 
cf)K + component, the remaining being pure phase space; the decay proceeds 
through the maximum estimated amount of fo{980)K + component, the re- 
maining being pure phase space. The estimated fractions, shown in Table 2, 
have been obtained by fits to the 3K invariant mass plots requiring that the 
K + K~ invariant mass lie within 2a of the nominal <fi mass for the <j)K + decay 
and between two kaon mass threshold and 1.05 GeV/c 2 for the /o(980)i^ + 
decay. These estimates are crude and represent conservative upper limits for 
the purpose of estimating systematic errors and are not meant to be measure- 
ments.0 Under these assumptions, the contribution to the total systematics 





D+ 


Df 


(j)K + 


12.4% 


18.75% 


fo(980)K+ 


44.5% 


72% 



Table 2 

Estimated fraction of (j)K + and /o(980)i"C + components for D + and Df decays. 

on the branching ratio measurement is ±0.10 x 10~ 4 for D + and X 10~ 3 
for D+. 

The last source of systematic error we studied is that due to fitting procedure. 
We calculated our branching ratios for various fit conditions, such as chang- 
ing the parametrization of the background shapes, rebinning the histograms, 
including in the D + fit the K~K + tc + reflection peaks and varying the fixed 
Df mass value by lcr of the quoted error [14]. Since all these results are a 
priori likely we used the resulting sample variance to estimate the associated 

2 We consider these as conservative upper limits since we do not account for the 
contribution of other components below the <f> and, when quoting the /o(980)-ftT + 
fraction, we do not simultaneously account for the (p. 
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systematics. We obtain a systematic contribution of ±0.19x10 4 for the D + 
decay mode and toll x 10~ 3 for the Df. 

In conclusion, summing in quadrature the different systematic errors we obtain 
our final results: 

BR (D+ -> K~K+K+) I (£>+ -> A-7T+7T+) = (9.49 ± 2.17 ± 0.22) x 10~ 4 
and 

5i? (D+ -> K-K+K+) I (D+ -> K-K+1T+) = (8.95 ± 2.12 ±|;§f) x 10~ 3 



4 Conclusions 

Our _D + measurement is consistent with the E687 upper limit [3] and consti- 
tutes the first clear evidence for this DCS decay. Our data indicate that only 
a minor fraction, if any, of the decay proceeds through the <pK + channel. This 
could suggest that the decay proceeds mainly through resonances that can 
couple to both irir and KK, such as the f resonance series, as expected from 
a naive spectator picture. However, more statistics would be needed to make 
quantitative statements through a Dalitz analysis. 

Our Df measurement is consistent with the E687 upper limit [3] and repre- 
sents the first observation of the 3K mode. It constitutes the second Cabibbo 
suppressed decay of the D+ measured. For Cabibbo suppressed decays other 
than Df — >• K + it~it + [12], only upper limits exist [13]. 
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